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ABSTRACT 


A  study  of  some  empirical  and  theoretical  functions  for  the 
effect  of  temperature  on  liquid  viscosity  indicated  that  the  viscosities 
of  the  n-alkanes  are  well  represented  hy  the  equation  log  (7? /d)  = 
log  (a/t)  +  b/(T-To),  where  d  is  the  density  and  A,  B  and  T©'  are  smooth 
functions  of  chain  length.  The  effect  of  increasing  static  pressure  is 
to  increase  A  and  To,  while  B  is  only  slightly  affected.  A  relaxation 
theory  developed  by  J.  H.  Gibbs  relates  B  to  the  barrier  to  rotation 
about  chain  bonds.  The  T©  (n)  values  are  predicted  for  n-alkanes  of  more 
than  n=6  carbon  atoms  by  the  Gibbs -DiMarzio  theory  with  a  flex  energy  of 
cal./mole,  in  agreement  with  estimates  by  other  methods.  A 
diagrammatic  method  for  evaluating  steric  effects  in  polymer  chains  was 
developed,  based  on  the  Pitzer  model  for  alkanes.  Some  excluded  volume 
effects  in  rings  and  short  chains  of  the  -CX2-  type  are  elucidated  in 
this  way.  An  empirical  approach  to  estimation  of  flex  energy  and  hence 
To  is  suggested.  A  statistical -mechanical  theory  of  the  glass  transition 
was  developed,  based  on  the  concept  that  the  observable  communal  entropy 
is  time  dependent.  At  Tg  a  contribution  to  the  communal  entropy  due  to 
backbone  motions  disappears,  giving  rise  to  a  time -dependent  glass 
temperature.  Glass  temperatures  are  theoretically  predicted  with  good 
accuracy  fran  viscosity  data  obtained  far  above  Tg.  The  Eyring 
transition  state  theory  for  relaxation  processes  was  modified  to  conform 
with  observed  polymer  behavior.  The  activation  free  energy,  enthalpy 
and  entropy  at  the  glass  temperature  were  evaluated  for  several  polymers. 
Intermolecular  interactions  in  liquid  alkanes  were  studied,  and  an  additive 
property  Q  =  (BpV)®*72^  was  found,  where  Bp  is  the  Antoine  vapor  pressure 
constant,  and  V  the  molar  voltane.  Q  is  dependent  on  the  number  of  C  atcms, 
the  number  of  side  branches,  and  the  number  of  pairs  of  vicinal  branches. 
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1.  Mechanism  of  Relaxation  in  Polymers 


The  temperature  dependence  of  the  viscosities  of  even  simple 
liquids  such  as  the  n-alkanes  cannot  be  adeqiiately  represented  by  the 
Arrhenixis  equation,  since  it  is  generally  found  that  the  activation 
energy  for  flow  Increases  as  the  temperature  is  reduced.  This  effect  is 
more  pronounced  in  high  poljmiers,  and  is  most  spectacular  in  the  transi¬ 
tion  zone  Just  above  the  glass  temperature.  MUler^  has  pointed  out 
that  the  empirical  equation  proposed  independently  by  Fulcher^  and  by 
Tammann  and  Hesse^, 


log  71  =  log  A  +  b/(T-To), 


(1) 


is  of  particular  interest  not  only  because  of  the  precision  with  which  it 
fits  accurate  viscosity  data  on  most  liquids  over  a  wide  tenqserature  range, 
but  also  because  it  can  be  transformed  into  the  WtP^  equation. 


log  aip  =  -Ci(T-Ts)/(C2  +  T-Ts),  (2) 


which  is  widely  xused  to  represent  the  temperature  dependence  of  relaxation 
processes  in  polymers.  In  these  equations,  7]  is  the  shear  viscosity, 

®T  ~  ^/7?s>  T|s  is  the  viscosity  at  a  reference  temperature  Ts,  and  A,  B 
To,  Cl  and  Cg  axe  adjustable  parameters.  It  can  be  shown  that  Cg  =  Ts-T© 
and  C1C2  =  B. 

Equation  (l)  indicates  that  7? - ^cx3  at  a  temperatiixe  T©  aot 

necessarily  equal  to  zero.  The  statistical -mechanical  theory  of  glass 
formation  developed  by  Gibbs  and  DlMarzio^  predicts  Just  this  type  of 
behavior.  It  is  therefoie  not  unreasonable  to  identify  T©  with  the  point 
at  which  the  configurational  entropy  of  the  liquid  goes  to  zero.  This  is 
estimated®  to  be  ItO  to  50*  below  the  glass  temperature,  and  can  be 
calculated  by  the  relationships  given  by  Gibbs  and  DlMarzio,  if  the  "flex 
energy"  (see  Section  3)  ia  their  theory  is  known.  The  purpose  of  tM« 
portion  of  the  program  was  to  determine  T©  for  chains  for  which  the  flex 
energy  is  known  from  independent  measurements,  for  coo^rlson  with  the 
theoretical  prediction. 
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The  compilation  of  selected  values  of  the  properties  of 
hydrocarhons  published  hy  the  American  Petroleimi  Institute's  Research 
Project  is  considered  the  most  reliable  source  of  precise  viscosities 
on  n-alkanes  over  a  wide  temperature  range.  The  dcmplete  set  of  data  on 
each  n-alkane  from  methane  to  eicosane  vas  fitted  to  Equation  (l)  by 
machine  computation  using  a  non-linear  least  squares  iteration  procedure, 
to  obtain  the  three  adjustable  parameters  for  each  cosqpound.  Additional 
values  vere  determined-for  n-alkanes  of  28,  36  and  6k  C  atoms  from  data 
published  by  Doolittle”.  ~ 

It  vas  fOTind  empirically  that  Equation  (3)  gives  a  consistently 


log  (TJ/d)  =  log  A  +  B/(T-To)  (3) 


better  fit  to  the  data,  as  can  be  seen  by  comparison  of  the  residual  sum  of 
squares  for  each  compound  in  Table  I.  d  is  the  liquid  density.  The  use 
of  kinematic  Instead  of  absolute  viscosity  is  indicated  by  Eyrlng's9 
treatment  of  viscous  flov  as  a  rate  process  and  by  the  statistical  theory 
of  elasticity  of  chalns^'O. 

The  precision  of  the  estimates  of  all  three  parsmeteirs  in 
equation  (3)  is  rather  poor,  even  thou^  the  viscosity  predictions  vhlch 
they  provide  are  quite  good.  This  is  because  of  the  hlc^  correlations 
among  the  parameters,  resulting  perhaps  from  the  fact  that  the  measurements 
vere  made  above  the  melting  point,  while  T©  is  very  much  lower.  Since  the 
data  vere  smoothed  by  the  American  Petrolermi  Institute  workers  before 
tabulating,  no  meaningful  statistical  statements  can  be  made  concerning 
the  precision. 

Gibbs  and  Ddlferzlo^  predicted  that  the  (free  energy)  barrier 
restricting  flov  from  one  configurational  state  to  another  should  Increase 
with  decreasing  temperature  as  T©  is  approached  from  above,  in  agreement 
with  the  Fulcher  eqiiatlon,  but  this  idea  was  expressed  only  qualitatively. 
Professor  Gibbs  has  recently  developed  a  quantitative  relaxation  theory, 
vhlch  he  advanced  during  his  last  visit.  Briefly,  the  theory  is  as  follows: 

A  segment  of  polymer  vhlch  is  in  some  high  energy  configuration 
as  a  result  of  the  application  of  a  shear  stress  will  relax  to  a  configuration 
of  lower  energy  at  a  rate 


kj  =  A  exp  (-zU/kT) 


ik) 


2 


where  z  is  the  number  of  bonds  in  the  segment  and  U  is  the  average 
potential  barrier.  At  least  two  configurations  must  be  available  in 
order  for  relaxation  to  occur,  so  in  general,  there  will  be  some 
minimum  configurational  entropy  required  for  relaxation.  Denoting 
this  by  s*',  and  observing  that  entropy  is  an  extensive  quantity. 


s*  =  SczVn 


(5) 


where  Sc  is  the  molar  configuration  entropy,  z*  is  the  number  of  bonds 
in  the  segment  of  critical  size  for  relaxation,  and  N  is  Avogadro's  niaaber. 
The  average,  rate  of  relaxation  for  the  entire  system  containing  segssents 
of  all  sizes  is 


A  exp  (-2zU/lsT)dz 

TT' 


exp  (-zU/l5T)dz 


A  /  -z^J  ' 

—  exp  - 

2  \  ^  I 


(6) 


Hence,  from  Equation  (5) 


A  -Ns*U 

k  =  - exp  -  .  (7) 

2  KTSc 


Assuming  that  the  configurational  entropy  is  responsible  for  the 
increase  in  specific  heat  at  Tg,  we  can  write 

rT 

Sc(T)  =  St-Sto  =  d  In  T.  (8) 

If  it  is  assumed,  as  a  first  approximation,  that  ACp  is  a  constant,  then 


Sc(T)  =  ACp  In  (T/To) . 


(9) 
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A 


Therefore , 


k  =  —  exp 
2 


-Hb^ 


KrAcp  In  (T/Tq) 


(10) 


Setting  the  relaxation  time  T  equal  to  k”^  and  combining  the  constants 
In  (10)  ve  have 


T=  A '  exp 


B' 

T  In  (T/To) 


(11) 


Professor  Glhhs  found  that  Equation  (11)  reproduced  the  V/LF 
function  over  a  vide  range  to  vlthln  a  few  per  cent.  The  Mathematical 
Analysis  Group  was  requested  to  test  the  equation  for  goodness  of  fit  to 
the  n-alkane  data.  At  the  same  time,  an  empirical  modification  of  these 
equations  vas  tested,  in  which  a  factor  of  T"l  is  Included  in  the  pre- 
exponentlal  tern.  The  equations  actually  used.  In  addition  to  EqTiatlons 
(1)  and  (3),  are  listed  below. 


log 

7J  =  log  (a/t)  +  B/(T-To) 

(12) 

log 

(7J/4) 

=  log  (A/T)  +  B/(T-To) 

(13) 

log 

7J  =  log  A  +  B/  T  In  (t/To) 

(Ik) 

log 

(7!/4) 

=  log  A  +  B/  T  In  (t/Tq) 

(15) 

log 

7}  =  log  (a/t)  +  B/  T  In  (t/Tq) 

(16) 

log 

iri/i) 

=  log  (A/T)  +  B/  T  In  (T/Tq) 

(17) 

The  sums  of  squares  of  the  deviations  from  these  functions, 
for  the  "best"  choice  of  the  adjustable  parameters,  three  in  every  case, 
are  shown  in  Table  I.  For  l8  out  of  the  23  liquids.  Equation  (13) 
had  the  smallest  deviation  from  the  observed  viscosities.  The  values  of 
the  parameters  determined  In  this  way  are  shown  in  Tables  II  throu^  V 
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TABIE  II 


Estimates  of  Parameters  in  Equation  (l)  for  n-Alkanes 


Number  of 


Carbon  Atoms 

log  A* 

B  (“K) 

To  (“K) 

1 

-1.533 

45.31 

37.0 

2 

-1.545 

136.6 

8.8 

3 

-1.409 

157.0 

21.3 

4 

-2.023 

458.8 

-68.7 

5 

-1.434 

204.5 

40.0 

6 

-1.537 

263.4 

36.3 

7 

-1.520 

273.8 

51.6 

8 

-1.634 

345.4 

41.2 

9 

-1.609 

346.8 

56.2 

10 

-1.698 

412.6 

46.2 

11 

-1.698 

427.9 

52.6 

12 

-1.749 

477.3 

46.2 

13 

-1.763 

501.0 

47.7 

l4 

-1.800 

542.3 

42.9 

15 

-1.797 

556.2 

45.8 

l£ 

-1.818 

588.8 

41.9 

17 

-1.794 

585.3 

49.1 

IB 

-1.818 

619.4 

43.9 

19 

-1.796 

616.8 

49.8 

20 

-1.798 

631.4 

50.3 

28 

-1.801 

771.4 

31.4 

36 

-1.605 

708.2 

64.3 

64 

-1.230 

659.9 

107.1 

*  Por  viscosity  in  cp 
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TABLE  III 


Estimates  of  Parameters  in  Equation  (3)  for  n^-Alkanes 


Nunflaer  of 
Carbon  Atoms 


log  A* 

-1.019 

-1.076 

-I.O65 

-1.073 

-1.130 

-1.132 

-1.196 

-1.191 

-1.235 

-1.236 

-1.256 

-1.259 

-1.271 

-1.265 

-1.270 

-1.254 

-1.258 

-1.238 

-1.237 

-1.249 

-1.121 

-0.829 


B  (”K) 

30.90 

90.70 

128.2 
261.8 

154.7 

337.8 
202.6 

242.7 

252.1 
286.6 

300.1 

325.5 

339.9 

360.1 

369.3 

385.7 

387.9 

402.8 
402.2 

411.4 

524.8 

514.6 

522.9 


To  (*1C) 

44.9 

25.8 

28.4 
-12.9 

54.5 

58.9 

71.3 

69.6 

81.2 

79.4 

85.4 

84.9 

88.1 

88.2 

91.8 

91.7 

96.9 

96.4 
101.6 

103.1 

86.1 

105.9 

133.7 


*  For  viscosity  in  cs 


TABIE  IV 


Estimates 

of  Parameters 

in  Equation  (l3)  for 

n-AUcaxies 

Number  of 

arbon  Atoms 

Iti  A* 

B  (“K) 

To  ("K) 

1 

3.0615 

18.132 

62.652 

2 

3.5813 

59.822 

56.084 

3 

3.6442 

142.73 

44.995 

k 

3.6763 

199.77 

44.833 

5 

3.9639 

152.57 

33.711 

6 

4.0272 

176.57 

97.581 

7 

4.0l^25 

215.33 

10Jk91 

8 

4.0632 

245.34 

113.29 

9 

4.0726 

278.69 

119.40 

10 

4.0896 

308.06 

124.76 

11 

4.1094 

334.79 

129.27 

12 

4.1319 

360.87 

132.51 

13 

4.1549 

383.41 

135.76 

14 

4.1771 

405.65 

138.37 

15 

4.2100 

420.65 

142.00 

l6 

4.2384 

436.85 

144.49 

17 

4.2773 

1^6.33 

148.30 

18 

4.3066 

458.88 

150.80 

19 

4.3520 

462.59 

154.99 

20 

4.3741 

475.06 

157.07 

28 

4.4808 

628.37 

lJi6.66 

36 

4.7276 

652.62 

158.78 

6k 

5.3347 

734.61 

174.92 

*  For  viscosity  in  cs. 
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TABLE  V 


Estimates  of  Parameters  in  Equation  (17)  for  n-Alkanes 


Number  of 


Dn  Atoms 

In  A* 

B  (  “K) 

To  f*K) 

1 

3.1167 

21.206 

61.1 

2 

3.6812 

89.775 

49.1 

3 

3.8519 

244.85 

36.0 

k 

- 

- 

5 

4.1277 

231.49 

73.0 

6 

4.1j831 

274.97 

83.6 

7 

4.2263 

328.52 

91.0 

8 

4.2390 

393.70 

94.7 

9 

4.2718 

435.67 

101.8 

10 

4.2841 

500.71 

103.4 

11 

4.3154 

541.18 

107.5 

12 

4.3381 

597.00 

108.3 

13 

4.3673 

635.47 

110.7 

14 

4.3899 

683.40 

111.4 

15 

4.4265 

705.97 

114.8 

16 

4.4541 

741.64 

115.7 

17 

4.4978 

747.92 

120.0 

18 

4.5242 

777.46 

120.9 

19 

4.5720 

771.80 

126.0 

20 

4.5958 

794.58 

127.4 

28 

4.7228 

1199.2 

105. 0 

36 

4.9945 

1165.7 

121.2 

64 

5.6579 

1213.6 

142.8 

*  For  viscosity  in  cs. 
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for  four  of  these  functions.  Except  for  Equation  (1),  all  of  these 
equations  yield  a  set  of  Tq  values  which  increase  smoothly  with  the 
number  of  carbon  atoms.  The  estimates  of  To  for  any  one  liquid  vary 
somewhat  from  one  equation  to  another,  and  there  seems  to  be  no  way  of 
choosing  among  them  on  grounds  of  physical  reasonableness.  Since  we 
are  primarily  interested  in  the  values  of  To  for  long  chains,  an  effort 
was  made  to  locate  viscosity  data  on  higher  polymers.  Tung^  reported 
zero  shear  melt  viscosities  as  a  function  of  temperatiire  for  branched 
and  xmbranched  polyethylene  fractions,  which  should  be  suitable  for 
calculation  of  To*  However,  the  points  were  too  scattered  for  a 
reliable  estimate . 

Aggarwal,  et  al.^  reported  viscosities  at  different  temperatures 
for  low  density  polyethylene  fractions,  and  these  were  fitted  to  Equation 
(1) .  Two  of  the  sets  of  data  were  too  scattered  to  be  useful,  and  Tq 
for  the  other  three  ranged  from  <<0  to  >>220*K. 

Marker,  et  al.^3  gave  melt  viscosities  for  a  whole  polymer  of 
low  density  polyethylene  with  Mn  =  33 >330.  The  data  fit  Equation  (1) 
with  To  »=  33*K  and  B  «  2,642*K.  These  values  do  not  seem  to  fit  in 
with  the  resxalts  for  the  n-paraffin  series.  It  was  concluded  that  the 
available  viscosity  data  on  polyethylehe  are  not  sufficiently  accurate 
for  our  purposes,  possibly  due  to  shear  degradation  or  oxidative 
degradation . 

,  While  this  work  was  in  progress,  the  paper  by  Karapet 'yants  and 
Tan^^  appeared,  giving  the  constants  in  Equation  (1)  for  35  hydrocarbons . 
However,  these  workers  foimd  it  necessaiy  to  use  a  separate  set  of 
constants  for  the  upper  (>0.85  cps)  and  lower  (<0.85)  ranges  of  values, 
making  a  total  of  six  adjustable  parameters  for  each  liquid.  Fxirthermore, 
many  of  their  Tq  values  are  negative,  in  contrast  to  those  in  Tables  II 
through  V  all  of  which  (except  for  butane)  are  positive.  Other  workers 
have  also  reported  negative  values  of  Tq.  This  seems  to  result  from  the 
use  of  viscosity  data  at  temperatures  far  above  To»  and  possibly  from 
Improper  weighting  of  data. 

The  viscosity  data  of  Fox  and  Flory^^  on  polystyrene  fractions 
can  be  tised  to  calculate  log  A,  B  and  T©  as  a  function  of  degree  of 
polymerization,  P.  The  data  had  already  been  fitted  to  an  equation  of 
the  same  form  as  (l)  by  Williams^,  and  as  Miller^  has  pointed  out,  the 
constants  in  Equation  (l)  can  be  obtained  by  a  simple  transformation.  The 
resulting  values  of  B  emd  To  are  shown  plotted  against  P  in  Figure  1. 

Both  B  and  To  approach  limiting  values  above  P  **  200,  and  drop  off  rapidly 
below  that.  This  type  of  behavior  is  predicted  by  the  Gibbs -DiMarzio 
theory,  but  no  atteii5)t  was  made  to  fit  the  data  to  the  theory,  since  Gibbs 
and  DiMeo'zio  had  already  obtained5  an  excellent  fit  to  the  glass  ten^jerattires 
of  these  san5)les. 
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2.  Pressiire  Coefficient  of  the  Glass  Temperatiire 


When  the  teii5)erature  of  a  liquid  is  reduced,  or  its  hydrostatic 
pressure  increased,  the  liquid  forms  a  glass  (unless  crystallization 
intervenes)  at  a  characteristic  temperature  which  increases  ^d.th  pressure. 

The  location  of  this  "glass  teii5)erature  curve"  on  the  phase  diagram  is  one 
object  of  this  study.  Glass  formation  is  preceded  by  a  marked  rise  in 
viscosity,  and  this  behavior  is  satisfactorily  represented  by  Equation 
(13),  for  exanqjle. 

An  atten5)t  was  made  during  this  program  to  develop  a  theory  for 
the  effect  of  pressure  on  viscosity,  to  complete  the  picture  phenomenologically. 
While  successful  to  a  degree,  the  adjustable  parameters  calculated  from 
published  data  were  complex  functions  of  temperature,  and  furthermore, 
significant  deviations  from  the  predicted  behavior  were  evident  on  plots 
of  log  7?  vs.  pressure.  On  inspection  of  families  of  isotherms  on  these 
plots,  iL  was  noticed  that  the  deviations  show  systematic  trends  with 
changing  temperature.  In  fact,  it  was  found  that  a  suitable  combination  of 
vertical  and  horizontal  shifts  results  in  almost  perfect  superposition  of 
all  the  isotherms.  As  an  example,  a  master  curve  is  shown  in  Figure  1  for 
1,  2,  4,  5»  6,  "Jf  8,  13,  l4,  15,  iS-dodecahydrochrysene,  in  which  all  the 

reported  datai7  for  this  liquid  are  s\iperposed  on  the  (extended)  135® 

Isotherm.  Two  additional  points  obtained  at  37.8"  are  not  shown,  since 
they  would  require  compression  of  the  viscosity  scale  to  accommodate  two 
more  logarithmic  decades,  although  they  also  superpose  nicely. 

The  pressTire  scale  for  each  isotherm  was  shifted  horizontally  by 
an  amount  pj  (atm. )  and  the  viscosity  scale  was  shifted  vertically  by  an 
amount  log  v^  to  produce  the  master  curve  in  Figure  2.  These  shift  factors 
were  found  to  be  well  represented  by  the  following  functions  of 
temperature: 


PT  =  42.1  (135-t),  (iB) 

log  vp  =  3.175  (IOOO/t-2.45),  (19) 


where  t  is  the  temperature  in  degrees  Centigrade  sind  T  is  the  absolute 
temperature.  The  coefficient  of  42.1  atm.  deg."l  in  Equation  (18)  has  the 
dimensions  of  a  thermal  pressure,  and  corresponds  to  a  cohesive  energy  of 
304  cal.  cc  at  25®.  The  activation  energy  from  Equation  (I9)  is  l4.5 
kcal.  mole”^. 
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In  order  to  determine  Tq  as  a  function  of  pressure,  it  is 
necessary  to  extrapolate  the  master  curve  in  Figure  2  to  log  v>j;  =00. 
Designating  the  corresponding  pressure  on  the  master  curve  as  P*(t), 
then  the  external  pressure  at  Tq  is  given  hy  Po*=P*(’t)-PTj  and  the  "To>  Pq 
curve"  for  dodecahydrochrysene  is  predicted  as 


Po  •=  P^(l35)-42.l(l35-t). 


(20) 


The  pressTire  coefficient  of  Tq  from  Equation  (20)  is  1/U2.1  =  0.0237 
deg.  atm.  "1,  which  is  at  the  upper  end  of  the  range  (0.004  to  0.024)  of 
values  for  ten  svihstances  tabulated  hy  O'Reillyi®  from  direct  observations 
of  glass  transitions. 

Of  the  fourteen  liquids  studied  by  I<nd.tz,  et  al.^’i^,  sufficient 
data  were  reported  on  seven  to  give  a  critical  test  of  the  general  validity 
of  the  tei35)erat-ure-Tpressure  superposition  principle.  It  was  thought 
desirable  to  vrork  ■vrf.th  kinematic  viscosities,  7?  /d,  since  it  had  been 
found  earlier  that  these  obey  the  Fulcher  equation  more  closely  than 
absolute  viscosities.  The  densities  of  some  of  the  liquids  were  given 
in  Reference  17,  but  had  to  be  interpolated  from  values  tabulated  at 
different  intervals  in  a  previous  publication^^  for  the  others.  The 
Mathematical  Analysis  Group  was  requested  to  perform  the  interpolations 
and  to  calculate  log  (77/d)  and  log  log  (7|/d).  R.  M.  King  provided  the 
author  with  values  of  d,  77 /d,  log  (77/^)  log  (77/d)  at  every 

tenperature  and  pressure  fOT  which  7|  'was  reported,  for  •the  seven  liquids 
chosen. 

It  was  found  that  the  log  (77/d)  -vs.  P  isotherms  for  perhydro 
chrysene  could  be  accurately  sijperposea  by  a  suitable  combination  of  -vertical 
and  horizontal  shift  factors,  Just  as  for  dodecahydrochrysene.  p^  was  found 
to  be  linear  in  tenperature,  and  log  -vj  ims  linear  in  T“l.  The  activation 
energies  and  pressure  coefficients  for  the  two  liquids  are  compared  in 
Table  VI,  where  it  can  be  seen  that  the  saturated  compound  has  the  lower 
activation  energy  and  higher  pressure  coefficient,  both  indicative  of  a 
lower  cohesive  energy. 
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TABLE  VI 


Effect  of  Teffi5)erature  on  Shift  Factors 


DodecahycLrochrysene 

Perhydrochrysene 


Activation  Energy 
(kcal.  mole-1) 

14.5 

9.3 


Pressure  Coefficient 
(deg.  atm. "1) 

0.0237 

0.0263 


The  precision  vith  which  the  viscosities  of  these . two  liqpiids 
fall  on  the  master  curve;  the  smooth  trends  in  shift  factors  with  changing 
tenroerature;  and  the  physically  reasonable  values  of  the  activation  energies 
and  pressure  coefficients,  all  lend  credence  to  the  principle  of  ten^jerature- 
pressure  superposition.  The  situation  is  actually  more  complex  than  this, 
however,  since  the  isotherms  for  the  other  five  liquids  studied  do  not 
superpose  on  a  plot  of  log  {TJ  /d)  vs.  P.  It  was  found  that  isotherms  for 
these  liquids  cotild  he  superposed  if  plotted  as  log  log  (7^/d)  vs.  P,  hut  it 
would  he  preferable  to  have  a  single  form  for  the  master  cturve,  especially 
if  a  good  analytical  representation  were  known. 

Since  Equation  (3)  is  known  to  fit  the  data  at  constant  pressure 
over  a  wide  ten^erature  range  on  all  the  liquids  for  which  it  has  been 
tested,  it  was  decided  to  construct  master  isobars,  which  could  then  be 
fitted  to  Equation  (3).  A  number  of  ftmctions  were  tried  before  it  was 
found  that  a  siniple  plot  of  log  7?/d  vs.  tes^jerature  forms  a  family  of 
isobars  which  exactly  superpose  when  shifted  by  appropriate  horizontal  and 
vertical  factors.  The  master  curve  is  of  the  form 


log  7|  /d  -  log  Vp  =  f(t-ep) 


(21) 


The  vertical  shift  factor  Vp  represents  the  effect  of  pressure  on  the 
parameter  A  in  Equation  (3)j  and  the  horizontal  shift  factor  0p  the  effect 
on  To.  The  fact  that  the  curves  for  all  pressures  superpose  means  that 
B  is  independent  of  pressure.  B  is  equivalent  to  C1C2  in  the  WEF  equation, 
and  it  was  originally  proposed20  that  this  is  a  universal  constant  with  a 
value  of  900°  for  all  polymers. 
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Master  curves  of  log  7|/d  vs,  t  reduced  to  one  atmosphere  were 
constructed  for  perhydrochrysene,  9"(2”cyclohe3{yethyl)heptadecane  and 
9“(2-phenylethyl)heptadecane  from  the  data  of  Lowitz,  et  al.l?  The 
precision  with  which  the  points  can  be  superposed  is  illustrated  for 
perhydrochrysene  in  Figure  3.  The  shift  factors  log  vp  and  6p  were  found 
to  be  smoothly  increasing  functions  of  pressure,  both  slightly  concave 
toward  the  pressure  axis.  Since 


To(P)  =  To(l)  +  ftp. 


(22) 


the  plot  of  8p  against  P  represents  a  phase  boundary  for  the  glass -liquid 
transition. 

The  Ehrenfest  condition  for  a  second-order  transition  is 


To“^dTo/aP  =  vA(X/ACp, 


(23) 


where  ACX  is  the  difference  in  thermal  coefficient  of  expansion  below 
and  above  the  transition  point,  and  the  difference  in  specific  heat. 
Since  both  thermal  expansion  euad  specific  heat  are  reduced  by  conpression, 
we  might  expect,  from  Equation  (23)  that  d  In  To/dP  might  be  independent  of 
pressure,  at  least  over  a  limited  range.  However,  it  was  found  that  the 
plots  of  log  To  vs.  P  were  also  concave  toward  the  pressure  axis. 

The  Mathematical  Analysis  Group  was  requested  to  determine  the 
best  least  squares  values  of  the  parameters  in  Equation  (3)  under  the 
restriction  that  B  be  independent  of  pressure,  for  the  four  liquids  for 
which  the  most  data  are  given  in  Reference  17.  A  computer  program  was 
written  for  the  simultaneous  estimation  of  the  least  squares  values  of  B, 
log  A  and  Tq  for  the  viscosities  reported  at  eleven  different  pressures,  a 
total  of  23  parameters  for  each  conpound.  At  the  same  time,  least -squares 
estimates  were  obtained  for  an  unconstrained  model  in  which  B  was  permitted 
to  vary  with  presstire,  giving  a  total  of  33  parameters  for  each  conpound. 
The  parameters  and  the  sums  of  squares  of  the  deviations  in  log  (7)/d)  for 
both  models  are  given  in  Tables  VII  to  X  for  the  four  liquids  studied.  The 
curve  drawn  in  Figure  3  represents  the  values  predicted  from  the  least- 
squares  parameters  in  Table  IX  with  B  constant.  Figure  4  is  a  plot  of  P  vs. 
Tq  for  all  four  conpounds. 
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TABIE  VII 

Effect  of  Pressure  on  Viscosity  of  9"(2*cyclohexylethyl)heptadeoaiie 
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TABUS  VIII 
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Effect  of  Presstire  on  Viscosity  of  Dodecahydrochrysene 
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Press\ire  ( atm. ) 


Effect  of  Pressure  on  T© 
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The  residual  sm  of  squares  is  lower  for  the  33  parameter 
unconstrained  model,  as  expected.  A  statistical  test  of  the  variance 
ratios  indicates  that  the  difference  is  significant  at  the  99.5^  con¬ 
fidence  level,  and  hence  that  a  single  value  of  B,  independent  of 
pressure,  is  not  supported  hy  these  data.  However,  it  should  he 
pointed  out  that  the  data  were  reported  at  integral  pressures  hy 
interpolation  of  the  raw  data,  so  that  the  calculated  variance  does  not 
give  a  true  indication  of  error.  To  establish  the  validity  of  the 
model  in  which  B  is  independent  of  pressure,  a  statistical  test  would 
have  to  be  made  on  the  original,  unsmoothed  data. 

The  fact  that  the  parameter  B  in  Equation  ( 3 )  is  only  a  slowly 
varying  function  of  pressure,  but  different  for  each  liquid,  suggests 
that  it  is  responsive  primarily  to  intramolecular  effects  and  much  less 
to  intermolecular  interactions.  According  to  Gibbs’  relaxation  theoiy 
(see  Equation  (lO),  Section  l),  B  »  Hs^/kACp.  Since  N,  s*,  k  and 
are  all  constants,  the  small  variation  in  B  with  pressure  must  be  due  to 
a  small  variation  in  U.  It  seems  reasonable  to  identify  H  with  the 
average  intramolecular  barrier  to  rotation  about  single  chain  bonds, 
which  should  be  relatively  insensitive  to  intermolecular  separation. 
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3.  Test  of  the  Gibbs -DlMarzlo  Theory 


The  purpose  of  this  section  is  to  test  the  Gibbs-DiMarzio 
theory^  against  the  To(n)  values  for  the  n-alkanes,  obtained  by  applying 
Equation  (3)  to  the  reported  viscosities,  and  to  estimate  the  flex 
energy  of  a  polymethylene  chain  by  this  means. 

The  Gibbs-DiMarzio  theoiy  considers  a  quasi -crystalline  lattice 
model  with  vacant  sites  (free  volume).  Every  bond  of  each  chain  is 
allowed  z*-l  possible  conformations,  where  z*  is  the  valence  of  the  chain 
atoms.  One  conformer  is  assumed  to  be  of  lower  energy  than  the  others, 
the  average  energy  difference  designated  as  G. 

The  energy  differences  calculated5  fi-cm  glass  temperat\ires  for 
several  polymers  have  values  which  are  reasonable  for  rotational  energy 
levels,  but  have  not  been  verified  by  independent  means.  The  corresponding 
values  of  C  derived  from  dilute  solution  viscosity  measurements  and  from 
thermoelastic  measurements  are  for  the  most  part  questionable,  as  pointed 
out  in  recent  conpilations  of  published  data21>22.  Only  in  the  case  of 
polyethylene  do  the  results  seem  to  be  in  ccoplete  accord  with  expectation, 
the  energy  difference  by  both  solution  cmd  bulk  measurements  being  about 
500  cal./mole23.  Furthermore,  this  is  the  e:]pected  magnitude  for  a  -CHs- 
chain  from  spectroscopic  measurements  on  butane,  pentane  and  hexane24. 

The  Gibbs-DiMarzio  theory  was  applied  to  values  of  To  calculated 
by  fitting  Equation  (3)  to  viscosity-tenperature  data  on  the  lower  n-alkanes, 
to  determine  whether  the  trend  with  increasing  chain  length  is  correctly 
predicted  with  a  flex  energy  G  of  approx±nately  500  cal. /mole. 

The  configurational  entrcpy  of  a  bulk  polymer  is  calculated  by 
the  quasicrystalline  lattice  method  developed  by  Huggins^?  and  Flory26 
for  the  entropy  of  mixing  polymer  with  a  low-molecular-weight  solvent. 

Solvent  sites  are  left  vacant  in  the  case  of  bulk  polymer,  corresponding 
to  a  free  volume  fraction.  The  mixing  is  considered  to  take  place  in  two 
stages:  (l)  disordering  of  the  perfectly  ordered  polymer,  and  (2)  mixing 
disordered  polymer  with  solvent  (vacancies).  It  is  necessary  to  take 
account  of  the  fact  that  there  is  an  energy  change  in  going  from  a  state 
of  perfect  order  to  one  of  random  arrangement.  There  is  no  way  to  estimate 
the  magnitude  of  the  change  in  intermolecular  energy,  so  the  energy  chaxige 
is  attributed  entirely  to  the  Intramolecular  effect  of  rotations  of  single 
bonds  into  higher  energy  conformations.  If  the  shape  of  the  chain  is 
sufficiently  synmetriccU.  when  in  its  lowest  energy  level,  then  crystallization 
is  predictedS?  to  occur  below  a  critical  value  of  f,  the  fraction  of  bonds 
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rotated  oat  of  their  ground  states.  If  not^  the  ziaaiber  of  configarations 
a-vallahle  to  the  system  irill  continue  to  decrease  as  the  tes^exature  is 
reduced  until  it  reaches  such  a  lov  value  that  farther  changes  occur  at  a 
negligible  rate  because  of  the  energy  barriers  separating  them.  The  system 
then  remains  frozen  in  one  configuration  at  lover  ten^ratures. 

The  relations  given  in  Reference  (5)  yield  the  foUoving 
expression  for  the  entropy  of  disorder  per  atom  vhen  chains  of  n  lattice* 
slte-occvqplera  of  valence  z*  are  placed  on  a  lattice  of  coordination  Timtiber 
z: 


S_  z-2  In  V  1  zv 

*  - - - -  +  —  (z-2  - - )  In 

R  2  1-v  2  1-v 


(z-2)n  +  2 


(1-v)  +  1 


znv 


+  n“^  In 


[(z-2)n  +  2]  (z-1)  n-3  , 

- +  - < 


ln[l  +  (z*-2)  g] 


(z*-2)  g 

+  - - 

1  +  (z*-2)  g 


(24) 


Here^  v  is  the  fraction  of  vacant  lattice  sites  (free  volume)^  g  »  exp 
(-CA^)  and  S  is  the  difference  between  an  average  energy  of  "flexed" 
bonds  and  the  energy  of  the  preferred  conformation. 

Equation  (24)  has  a  nuaiber  of  interesting  properties  which  should 
be  pointed  out.  Rhen  ^z*»4  and  S-  is  set  equal  to  zero  to  solve  for 

negative  values  are  obtained  for  ns^  or  less^  meaning  that  random 
p^^v^Tle  ia  always  possible  no  matter  how  stiff  the  chains,  provided  they 
are  shog^ter  than  six  units.  ,  Since  we  have  To  data  for  chains  of  as  few 
as  five  atoms,  we  assume  that  t;te  lattice-site  occxpler  is  a  single 
metbylene  group,  and  seek  to  reduce  the  entropy  by  varying  some  of  the 
parameters  in  Epiatlon  (24).  Adjustment  of  the  free  volume  v  over  a 
reasonable  range  of  values  does  little  to  improve  the  fit  to  the  data. 

Hence  it  is  expedient  to  neglect  free  volume  in  calculating  the  flex 
energy  C  . 
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Gibbs  and  DiHarzlo  assumed  that  the  conformational  energy  of 
each  bond  is  independent  of  the  conformations  of  nei^boring  bonds.  In 
actuality^  certain  ccnibinatlons  of  torsional  angles  result  in  large  steric 
effects.  In  the  case  of  polymethylene  chain^  in  particular^  a  gauche  bond 
foUoved  imnediately  by  a  gauche  z^^tation  in  the  opposite  direction  (G6* ) 
is  virtually  excluded  by  reason  of  its  high  energy.  Mathematical  techniques 
for  treating  the  pr^lem  of  pairwise  correlations  of  rotation  have  be«Ei 
developed  by  Llfsoi^  and  applied  to  the  calculation  of  random  coil  dimensions 
by  Hagai  and  Ishikawa^  and  by  Hoeve30^  and  to  the  entropy  of  melting  by 
Starkweather  and  Boyd3i.  The  decrease  in  entrc^  due  to  the  exclusion  of 
GG*  pairs  (pentane  effe^)  is  ^eatest  when  *»  0^  and  is  negli^ble^^31 
for  CA  ^  5  deg."^  mole”!.  For  n->oo.  Expiation  (2U)  gives  CAo  *  3.81 
(for  v«0  and  z«z**»4),  and  GAo  increases  as  n  ^creases.  The  formula  given 
Tty  Taylor 32  for  the  fraction  of  "forbidden"  confonaations  yields  0,05  R 
for  the  entropy  decrease  in  pentane^  and  this  is  too  small  to  resolve  the 
difficulty  with  the  calculation  of  To.  Hence^  modification  of  the  Gibbs- 
DiMarzio  theory  to  include  pairwise  correlations  does  not  seen  warranted. 

Setting  Sc>bO  in  Equation  (24)  and  letting  we  arrive  at  the 
final  expression  to  be  cosqpared  with  the  set  of  (To»  u)  values  for 
kinematic  viscosities  listed  in  Table  III: 


Since  Equation  (25)  predicts  a  negative  To  for  im>5^  it  is  not 
surprising  to  find  a  large  error  in  the  prediction  for  iv^.  Therefore^ 
only  the  values  of  To  for  np*7  through  20  were  used  in  the  cooqnztatlon  of 
the  least -squares  value  for  which  was  found  to  be  490.8  cal./aole.  The 
calculated  values  of  To  are  shown  in  Table  XI  for  this  choice  of  C .  The 
agreement  with  To  calculated  from  viscosities  is  considered  excellent,  tlw  , 
flex  energy  of  490.8  cal. /mole  agrees  with  that  estimated  by  other  mcthods^3,24^ 
and  the  limiting  value  of  To  ■■  1^.8*K  as  n  ecpproaches  oo  seams  reasonable 
for  linear  polyethylene. 
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TABIE  XI 

Gibbs -DiHarzlo  Equation  Using  Bata  From  Heptane  to  Elcosane 

G  »  490.8 


n 


Catleulated  Tp 


(Observed  Tp 


7? 


-  Calculated  T 


100.1 

3.0 

99.1 

2.5 

97.9 

-1.5 

96.6 

0.3 

95.2 

-3.5 

93.6 

-1.8 

91.8 

-3.6 

89.7 

-1.6 

87.3 

-2.4 

84.5 

0.9 

81.0 

-1.6 

76.7 

4.5 

71.0 

-1.4 

62.6 

8.7 

as  n - >00*  128,8 
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Glass  faroatlon  has  of  course  not  hean  observed  In  the  n-aUcanes 
because  of  the  hi£^  rates  of  crystallization  in  these  coopounds^  but 
Cohen  and  Tuniball33  point  out  that  there  is  reason  to  believe  that  any 
liquid  vould  fom  a  glaiss  if  cooled  rapidly  enough.  Naturally^  there  is 
no  available  data  on  the  properties  of  these  glasses^  so  that  ve  cannot 
independently  evaluate  the  volume  of  a  lattice  slte^  or  the  van  der  Waals 
energy^  which  are  parameters  in  the  Gibbs‘*DlIfeLrzio  theory.  Veglect  of 
free  volume  is  equivalent  to  assuming  a  van  der  Waals  energy  which  is 
large  conpared  to  the  flex  energy.  Ihis  is  probably  justified  in  the 
present  case. 

The  use  of  a  lattice  model  for  liquids  is  in  many  respects  a 
poor  approximation^  although  the  approximation  becomes  more  appealing 
below  the  melting  pointy  and  especially  as  the  glass  tenperature  is 
approached.  The  fact  that  a  reasonable  prediction  is  made  for  heptane 
and  all  the  higher  honologues  is  a  remarkable  confirmation  of  the  Glbbs** 
DiMarzio  theory  of  glass  formation. 
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4.  Tgg-faiwfl.tlon  of  Flex  Energy 


The  relative  eixergies  of  different  rotational  confomations 
(confonoers)  are  required  in  order  to  apply  the  Gibbs -DiMarsdo  theory 
to  prediction  of  glass  teopexatures.  This  problem  also  arises  in  the 
interpretation  of  n.m.r.  spectra,  infrared  spectra,  chain  dimensions  of 
polymers  in  solution,  and  stress Hjemperature  coefficienbs  for  cross- 
linked  rubbers.  Most  investigations  have  assroned  that  the  main  contributing 
factor  to  the  rotational  potential  is  the  interaction  between  non-bonded 
atoms  adjacent  to  the  bond  about  which  rotation  takes  place.  Ifoson  and  ' 
Kreevpy^,  for  example,  calculated  the  rotational  potentials  for  substi¬ 
tuted  ethanes  from  the  interactions  betweei  substituents.  The  interaction 
potentials  were  approximated  by  those  found  for  rare  gases  having  similar 
electronic  structrrre  to  the  substituent.  For  example,  the  interaction 
between  F  atoms  was  replaced  by  the  potential  found  from  the  second  virial 
coefficient  of  neon.  A  correction  was  added  for  the  interactions  between 
bond  dipoles.  These  same  potential  functions  were  later  used  by  Idquori, 
et  al.33,  in  calculating  the  configurational  energies  of  isolated  polymer 
chains. 


In  the  calculations  discussed  above,  the  known  bond  lengths  and 
angles  were  used.  These  are  very  difficult  to  predict  for  aiy  given  coapound, 
yet  the  non-bonded  Interactions  are  quite  sensitive  to  their  exact  values. 
Even  when  the  known  lengths  and  angles  are  used,  the  calculated  energies  are 
in  only  rough  agreement  with  esqperlment.  It  has  been  suggestedS^  that  in 
addition  to  the  non-bonded  interactions,  each  bond  has  an  "intrinsic" 
resiistance  to  torsion.  The  contribution  of  this  effect  to  rotatloiaal 
energies  is  not  known. 

In  view  of  the  rmcertainties  in  the  origin  of  the  energy 
differences  between  conformers,  it  is  not  surprising  that  the  calculated 
values  are  very  unreliable.  Nevertheless,  even  a  crude  estimate  of 
conformational  preferences  would  be  of  assistance  in  estimating  flex 
energies,  or  at  least  in  ranking  structures  in  order  of  increasing  stiffness. 
The  following  simple  model  for  steric  effects  in  alkanes  was  introduced  by 
Pitzer37  in  1940.  The  carbon-carbon  bonds  are  assumed  to  be  of  equal 
length,  to  meet  at  tetrahedral  angles,  and  to  be  in  staggered  conformations. 
Certain  configurations  can  be  excluded  from  consideration  because  of  atcm 
overlap.  Certain  others  result  in  close  approach  of  one  or  more  pairs  of  H 
atoms,  and  these  are  assigned  an  energy  a  per  pair.  With  a  «  800  cal. /mole, 
the  model  was  foimd  to  account  satisfactorily  for  entropies  determined 
spectroscopically. 
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Pitzer  used  Fischer-Hirschfelder  models  to  determine  the  manber 
of  interactions,  but  models  can  be  very  umrieldy  when  the  number  of  atoms 
to  be  considered  is  large.  During  this  program  a  scheme  was  devised  for 
diagraomatically  representing  spatial  configurations  without  atom  overlap, 
so  that  the  number  of  non-bonded  nearest  neighbors  can  be  determined  for 
the  various  allowed  configurations,  and  hopefully  correlated  with  known 
values  of  chain  stiffness.  The  atoms  are  located  on  a  tetrahedral  lattice, 
shown  in  projection,  with  the  relative  levels  of  bonded  atoms  depicted  by 
the  "flying  wedge"  syaibol.  The  bonds  normal  to  the  plaxie  of  the  projection, 
so  that  the  atoms  bonded  are  superimposed,  appear  as  in  a  Newman  projection. 
Only  staggered  confoxmatlons  are  possible  on  a  diamond  lattice. 

Schematic  representations  of  some  polymethylene  chains  are  shown 
in  Figure  3,  The  bonds  are  designated  trans(T),  gauche  positive  (g),  or 
gauche  negative  (G* )  depending  on  the  orientation  of  the  preceding  azul 
following  bonds.  Tte  allH:rans  sequence  cont€dns  no  non-bonded  nearest 
neighbors.  Rotation  of  one  bond  by  120*  to  the  Muche  position  introduces 
one  pair  of  non-bonded  nearest  neighbor  H  atoms  (structure  I).  An  additional 
G  bond,  either  adjacent  to  the  first  (structure  H),  or  separated  ftrom  it 
by  one  T  (structure  IH),  introduces  an  additional  pair.  Hence  the  number 
of  such  pairs  of  H  atcais,  with  interaction  oiergy  a,  is  simply  equal  to  the 
nuidber  of  G  bonds.  This  is  a  justification  for  the  use  of  a  single  value 
of  the  "flex  energy"  for  all  of  the  n-alkaaes. 

The  pairs  of  Intenctlng  H  atoms  occtpy  adjacent  positions  on  an 
open,  slx-mesibered,  cyclohexane-chalr-type  ring,  which  can  appear  on  the 
diagram  In  only  two  ways:  on  edge  as  in  structure  I  or  flat  as  in  U  and 
III  (which  have  the  additional  ring  on  edge). 

Configurations  which  are  disallowed  because  of  atom  overlap  are 
easily  detected  if  the  overlapping  atoms  are  sQMrated  by  six  bonds,  since 
these  occupy  the  same  position  on  a  closed  six‘inaia>ered  ring.  More  distant 
atoms  may  overlap,  and  It  was  decided  to  study  the  conditions  under  which 
this  occurs.  It  was  pointed  out  by  Taylor 32  that  for  chains  of  the  type 
-CXs*,  where  X  is  a  monovalent  atom,  overlap  of  X  atosui  separated  by  six 
bonds  occurs  when  any  pair  of  successive  chain  bonds  have  the  conformation 
GG*.  Nagai  and  Ishikava29  determined  that  the  next  longest  disallowed 
sequence  for  this  type  of  chain  is  GGTGGTGG  or  its  mirror  image  (formed  by 
ireplaclng  every  G  with  G').  This  sequence  results  in  double  occupation  of 
one  site  on  a  twelve-memibered  zing,  as  shown  for  unde  cane  in  Figure 
structure  17.  In  general,  the  number  of  sequential  rotations  required  to 
produce  double  occtpatlon  of  a  ring  site  equals  the  ring  size  minus  four, 
for  -CXa-  type  chains. 
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Intramolecular  Interactions  In  -CXa**  Chains 


Overlap  to  close  a  fourteen-membered  ring  vas  found  by  the  use 
of  the  dlagraauB  to  be  possible  id.th  six  different  bond  sequences  not 
Including  GG*  or  GGTGGTGG.  An  additional  six  sequences  vhich  close  a 
8ixteen-«ena>ered  ring  were  also  found.  The  study  was  discontinued  at 
this  pointy  partly  because  of  the  publication  by  Sadth^^  of  the  enuaera- 
tlon  of  disallowed  sequences  for  rings  of  iq>  to  eighteen  aeBbers,  obtained 
with  a  coBQnzber.  This  reference  should  be  consulted  for  application  of 
the  nethod  to  the  problens  of  chain  dimensions  and  ring  stability. 

The  above  discussion  concerns  only  "one  atom"  chains  of  the 
type  -CXa-  where  X  is  a  single  atom.  "Two  atom"  chains  such  as  KJHbKIHX-, 
or  *<IH2**CX2'*  would  be  of  greater  interest,  but  involve  farther  cc^plications: 
first,  the  pairs  of  side-chain  atoms  separated  by  five  bonds  can  be  of  two 
types,  R— X  or  H— H,  and  second  the  sterlc  effects  resulting  frcm  bond 
rotation  are  probably  no  more  important  than  the  electrostatic  effect 
resulting  from  the  change  in  orientation  of  bond  dipoles.  In  the  case  of 
poly( vinyl  chloride)  (PVC)  the  electrostatic  interactions  have  been 
estimated  by  Fordhai^^  by  the  use  of  a  diamond  lattice  model,  with  the 
results  summarized  in  Table  XII.  The  eight  possible  configurations  of  the 
two  bonds  to  a  methylene  group  in  Isotactic  (i)  and  syndlotactlc  (s)  PVC 
resTilt  in  four  different  electrostatic  intearactlon  energies.  The  sterlc 
interaction  energies  from  Pltzer's  model  are  included  in  Table  XU.  The 
most  stable  configuration  for  s-FVC  is  predicted  to  be  TT,  in  agreement 
with  experiment,  the  next  highest  energy  configuration  is  T6,  and  the 
difference  in  energy  is  predicted  to  be  67O  +  a.  The  corresponding 
predictions  for  i-FVC  are  TG,  GG,  and  again  a  dlffeience  of  67O  a. 

The  value  of  To  for  high  molecular  weight  PVC  is  seen  in  Table  XU, 
Section  6  to  be  42*C.  As  shown  in  Section  3^  the  Gibbs -DiMarzio  theory 
predicts  that  €  »  3.81  To  for  long  chains.  Setting  this  equal  to  67O  +  a, 
a  value  of  To  “  33£*K  corresponds  to  a  »  5^0  cal. /mole,  which  is  close  to 
the  H— H  interaction  in  alkanes.  The  H— Cl  interaction  which  occurs  in 
alternate  TG  configurations  is  presumably  of  higher  energy  and  does  not 
contribute  appreciably  to  the  glass  transition. 

While  the  results  of  the  calculation  just  described  are  in 
satisfactory  agreement  with  experiment,  the  method  has  not  yet  been 
extended  to  other  polymers  because  of  additional  problems  which  arise. 

For  exanple,  the  most  stable  configuration  of  pol^vlnylidene  chloride) 
(PVCI2)  is  TGTG*^®.  This  sequence  is  determined  by  the  relative 
orientations  of  three  successive  methylene  chloride  grojqjs,  requiring 
evaluation  of  steric  and  electronic  energies  of  six  allowed  configurations 
on  a  diamond  lattice.  Even  longer  sequences  should  be  examined  in  the  case 
of  polymers  whose  configurations  have  not  been  established,  to  insure  that 
the  configuration  of  minimum  energy  will  be  found. 
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Electrostatic  axid  Sterlc  Interaction  Energies  in  FVC  (cal./iaole) 
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Another  prohlem  which  has  not  "been  solved  concerns  the  fact 
that  poisoners  with  two  bTiUsy  substituents  on  the  same  chain  atom,  such 
as  polyisobutylene,  have  no  aJJLowed  configurations  on  a  diamond  lattice. 
This  will  require  assignment  of  interaction  energies  to  pairs  of  atoms 
occupying  the  ssune  lattice  site,  to  determine  the  configuration  of 
lowest  energy. 

The  problem  of  energy  differences  between  rotational  isomers 
in  polymer  chains  has  been  studied  by  various  means,  all  dependent  on 
the  fact  that  the  mean  square  separation  of  the  chain  ends  is 

determined  by  the  energy  differences  and  hence  varies  with  teng)erature. 

In  recent  reviews21>22  of  published  values  of  d  ln<ro2>/dr,  it  was 
concluded  that  while  satisfactory  agreement  can  be  reached  using 
different  methods,  no  satisfactory  interpretation  of  the  results  exists. 

In  particular,  the  theory  developed  by  Ptitsyu?3  predicts  zero  or 
negative  values  of  d  ln<Cro%>/dT  for  all  vinyl  polymers,  whereas  this 
coefficient  is  actually  found  to  be  positive  in  most  cases. 

The  dilute  solution  values  of  d  Incfro^^/dT  are  open  to  some 
question,  due  to  specific  polymer-solvent  interactions.  Thermoelastic 
measurements  on  bulk  polymers  are,  of  course,  not  subject  to  this 
perturbation,  and  hence  are  probably  more  reliable.  The  flex  energies 
computed  from  these  measurements  shoxild  agree  with  the  Gibbs -DiMarzio 
estimates  of  flex  energy  from  glass  temperatures.  However,  it  should 
be  pointed  out  that  the  chain  dimensions  are  very  sensitive  to  correlations 
of  rotations  of  neighboring  bonds,  and  the  flex  energies  calculated 
depend  tpon  what  assumptions  are  made  concerning  bond  angles,  rotation 
angles  and  correlations  of  rotations. 

Table  XIII  shows  flex  energies  C  calculated  by  Ptitsyn23  from 
thermoelastic  measurements,  except  for  the  value  for  isptactlc  poly(methyl 
methacrylate),  which  was  calculated  by  Sakurada,  et  al.*^^  from  intrinsic 
viscosities  in  different  theta  solvents.  The  corresponding  flex  energies 
from  the  Gibbs -DiMarzio  theory  shown  in  Table  XHI  were  calculated  from 
the  relation 


€  »  3.81  To  (^)  (26) 


as  discussed  in  section  3.  The  Tq  values  for  polyisobutylene  and  natural 
rubber  were  taken  from  Perry*  s  book^,  using  the  relation  Tq  “  Tg-Ca,  where 
Ts  is  the  reference  temperature  (Tg  in  these  cases).  For  the  silicone. 
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To  vw  estiiaated  from  Ts  and  Cz  reported  by  BarloiA^^  et  al.  On  the 
assunption  that  Cz  is  independent  of  tacticity^  To  for  isotactic 
poly( methyl  methacrylate)  vas  calculated  from  the  glass  tenperature  of 
*♦8*  for  the  isotactic  polymer,  and  Cz  "  44"  reported  by  Saito,  et  al.43 
for  conventional  polymer.  The  value  of  To  for  polyethylene  is  taken 
from  section  3#  Table  XI. 

The  agreement  seen  in  Table  XIII  for  the  first  three  polymers 
is  encouraging,  but  for  polyisobutylene  the  two  methods  disagree  as  to 
ojrder  of  magnitude.  Moreover,  if  Ptitsyn*s23  equation  for  isotactic 
vinyl  polymers  is  used  to  calculate  flex  energies  from  the  data  coomiled 
by  Ciferri22  for  various  polyacrylates,  the  largest  G  found  is  97*^ 
cal. /mole  for  poly( ethyl  acrylate),  whereas  G  ■  970  from  Equation  (26) 
with  To  “  Tg-40  »  The  atactic  nature  of  the  polymers  probably 

accounts  for  the  smedl  effect  of  temperature  on  chain  dimensions  of  the 
polyacrylates.  That  is,  even  in  its  lowest  energy  coixfiguration  an 
atactic  polymer  will  be  highly  coiled  due  to  the  random  arrangement  of 
d  and  1  tinits,  which  differ  in  directional  properties.  The  higher 
energy  states  are  not  necessarily  smaller  in  <  ro%>  • 
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5.  The  Cell  Model  of  the  Glass  Transition 


A  simple  kiiietic  model  for  the  glass  tztmsltion  has  been 
formulated  using  the  conc^t  of  communal  entropy.  Communal  entropy  is 
a  contribution  to  the  total  entropy  of  a  system  vhich  exists  trhen  the 
molecules  or  other  elements  of  the  system  are  indistinguishable^.  The 
amount  of  communal  entropy  in  the  system  is  determined  by  the  partitioning 
of  the  system  into  cells  containing  indistinguishable  molecules^  l.e.^  the 
degree  to  vhich  the  molecules  in  the  system  are  Interchangeable.  If  there 

are  n  molecules  per  cell,  the  volume  per  cell  is  V/-q-,  and  the  nosiber  of 

permutations  of  n  molecules  in  all  N/n  cells  is  ;Just  (ni  )®/^.  Hirschf elder, 
Stevenson,  and  Eyring^5  felt  that  solids  acquire  essentially  tte  total 
communal  entropy  on  melting.  This  viev  vas  criticized  by  Rice^  and 
Gurney  and  Mott^7.  The  latter  authors  showed  that  a  liquid  possesses 
almost  2/3  of  the  full  communal  entropy  with  only  5  laolecules  per  cell. 

If  it  is  assumed  that  the  main  glass  transition  in  polymers  is 
due  to  motions  involving  s  backbone  atoms,  vheire  s  «  ^”10,  it  is  possible 
to  construct  a  simple  model  for  the  glass  traiusition  based  on  the  communal 
entropy  concept.  The  glass  transition  should,  occur  when  there  are  s  atoms 
per  cell. 


In  a  system  containing  N  polymer  molecules,  there  vUl  be  a  term 

in  the  partition  function  eqpal  to  l/C*^)!,  tthere  there  are  n  chain  atoms 

per  cell.  As  the  temperature  is  lowered  toward  this  term  will  remain 
unchanged  in  the  equilibrium  partition  function.  However,  measurements 
made  using  a  finite  time  scale  will  reflect  a  gradual  decrease  in  n.  At 
Tg,  n!»s,  and  there  will  no  lon^  be  an  entropy  contribution  from  the 

l/(-^)l  term.  Thus  by  producing  an  expression  for  ~  in  terms  of  T, 

it  should  be  possible  to  evaluate  Tg. 

If  we  partition  all  of  the  volume  into  cells  of  n/s  segments, 
the  number  of  segments  per  cell  is  given  by: 


n  vPn 

»— •  &  - 

s  M’s 


(27) 
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irhere 


V  ■  the  cell  voixnne 
N  *»  Avogadro's  nuuiber 
P  *=  the  polymer  density 
M*  *»  the  molecolar  ireight  per  hackhone  atom. 


The  cell  volume  v  is  obtained  as  the  volume  of  the  sphere 
vhose  radius  is  equal  to  the  root  mean  square  distance  travelled  by  each 
segment  in  time  t  In  self-diffusion. 


4 

V  ts  — - 

3 


TT  (6itetD)^^^ 


(28) 


Here  D  is  the  molecular  self -diffusion  coefficient  and  na  is  the  number 
of  chain  atoms  per  molecule.  Equations  (27)  and  (28)  are  combined: 


n  pil(4-7r)(6n£tD)3/2 

«—  a  - - « - - 

6  M's 


(29) 


This  expression  gives  the  ntonber  of  segments  per  cell  in  terms  of  D, 
Unfortunately,  there  are  very  few  data  reported  in  the  llterattrre  for 
self -diffusion  of  polymers.  Bueche^  has  derived  an  expression  relating 
D  to  the  viscosity  coefficient,  TJ  : 


pmPr 

36m7;^  ' 


(30) 
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where 


R  is  the  gas  constant,  8,3l4  x  107  ergs  mole”^  clfig”^, 

h^  is  the  tugaerturbed  ]nean<>square  end-to-end  distance 
in  cnu,  and 

M  is  the  pol^ioer  molecular  weight. 


Substituting  into  Equation  (29),  we  obtain 


JL  .  tRli%  \  3/2 

s  \  M-  1  i  3s  /  1  671  I 

The  viscosity  is  well  r^resented  by  the  Fulcher^  equation: 


(31) 


B 

log  7]  ss  log  A  +  - - -  ,  (32) 

T-To 


The  appropriate  substitution  is  made  to  reach  the  final  expression: 

P  \  /  IfTTN 

M’  /  \  3s 

The  condition  for  is  that  — —  b  l.  Thus,  at  T„ 

g  ^  '  g 

5  Ip]  2  HTTn  tRh^T 

0  *s  —  log  I  — —  I  +  — —  log  .  +  log  - -  log  A 

3  \  M’  /  3  3s  6 
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Equation  (3^)  is  a  transcendental  equation  and  therefore  a 
solution  cannot  he  obtained  hy  oz‘dlnary  algebraic  iaethods.  However^  a 
solution  is  obtained  rapidly  by  using  a  self-consistent  iterative  method. 
Since  Tg  «  0  (  300*K),  it  follows  that  !«!>>> log  Tg,  Therefore,  if  an 
i^roximate  v^ue  of  Tg  is  used  in  the  log  Tg  term,  the  equation  becomes 
linear  in  Tg.  This  equation  is  solved  for  and  then  the  solution  is 
used  as  the  estimate  in  the  log  Tg  texm.  The  process  is  repeated  until 
the  solution  and  estimate  agree  within  the  desired  tolerance.  Ordinarily 
no  more  than  two  iterations  are  required. 

Althou^  Equation  (34)  is  the  result  of  a  rather  simple 
deilvation,  it  may  be  seen  readily  that  it  is  in  qualitative  agreement  at 
least  with  current  theories  of  the  glass  transition.  The  equation  predicts 
that  Tg  will  increase  with  the  monomer  molecular  weight  and  decrease  with 
increasing  time  scale  of  the  experiment,  approaching  a  lower  limit  of  To 
as  t  approaches  infinity.  The  latter  prediction  agrees  with  the  Gibbs 
and  DiMarzlo^  theory  of  a  lower  limit  to  Tg  at  an  equilibrium  second-^rder 
transition  tcnperature,  T2.  The  inverse  relationship  between  Tg  and  the 
end-to-end  distance  h^  in  Equation  (34)  can  be  disturbing  imless  one 
considers  the  effect  of  the  viscosity  (see  Equation  31).  It  is  well-known 
that  in  the  high-molecular  weight  region5 


<h^o  »  am  (35) 


and 


(36) 


Thus  the  net  result  is  a  positive  dependence  of  Tg  on  M. 

A  test  of  the  theory  was  made  using  Miller's^  results  for  two 
lowiiralecular-wei^t  polystyrene  fractions,  for  the  choices  s*5>  t*»1000  sec., 
and  Pssl.O  g/cc.  The  results  are  given  below  in  Table  XIV. 


I4O 


TABIE  XIV 


Glass  Tes^>eratures  of  Polystyrene  Calculated  with  Equation  (34) 


Myls 

Tc,  %  Exptl. 

Tg,  “K,  Eq.  (34) 

Error,  “K 

66^ 

350 

353 

+3 

3041 

337 

336 

-1 

1675 

313 

319 

46 

Viscosity-tenperature  data  based  on 
fitted  to  an  equation  of  the  form  of  Equation 


a  mmiber  of  polymers  were 
(32), 


B 

log  /)  T  «  log  A'  +  - - ,  (12) 

T-To' 


in  order  to  obtain  activation  parameters.  (This  will  be  discussed  in  more 
detail  in  the  next  section. )  The  results  weire  used  to  calculate  glass 
temperatures  for  these  polymers  after  Equation  (34)  was  modified  to  conform 
with  Equation  (12): 

5  I  ^  \  2  tEh%g2 

0  »  —  log  — —  +  —  log  +  log  "  . -  -  log  A* 

3  \  M*  /  3  3s  6 


-  S'/iTg-To*) 


(37) 


Glass  tsa|;>eratTires  vere  calcalated  vLth  tslOOO,  sod  P=1  as  before. 
All  values  of  In  this  work  vere  taken  fjrom  Kurata  and  Stockmayer's^^ 
values^  vhlch  are  based  on  intrinsic  viscosities.  The  results  are  given 
in  Table  XV. 

In  general  it  aagr  be  seen  that  the  results  are  in  good  agreement 
vlth  experiment.  Tbds  is  especially  gratliying  in  view  of  the  fact  that 
the  lowest  measurement  teaperature  was  often  veil  above  the  experimental 
glass  teaperature.  In  the  ease  of  poly(decaffletfaylene  adipate)^  for 
instance,  the  glass  teaperature  is  predicted  vlth  an  error  of  6*K  trcm  a 
point  above  it. 

During  this  study  it  vas  found  that  virtually  all  recent 
publications  concerned  vlth  transitions  in  polyisobutylene  use  the 
Flory-Fooc  equation^5^ 


Tg,  *0  »  -63  -  6.9  X  IC^/Af 


(38) 


to  espress  the  i*elatlonshlp  between  Tg  and  akolecular  weight.  Equation  (38) 
is  in  fact  based  on  the  measurement  of  Tg  for  a  single  polymer  saaple^S, 
plus  certain  free-colnme  considerations.  The  few  other  values  of  Tg 
appearing  In  the  literature  do  not  conform  closely  to  this  equation^  there¬ 
fore^  only  one  experimental  value  of  Tg  Is  given  in  Table  ICV  for 
polyisobutylene . 

It  vas  demonstrated  In  this  program  that  the  use  of  kinematic 
viscosity  instead  of  absolute  viscosity  in  the  Fulcher  equation  results 
in  a  better  fit  of  the  experimental  data: 


log  TJ  m  log  P  +  log  A  + 


B 

T-To 


(39) 


This  of  course  results  in  changes  In  each  of  the  three  parameters  A,  B, 
and  Tq.  Substitution  of  Equation  (39)  luto  Equation  (31)  leads  to  the 
following  exqoresslon  for  Tgt 


Tg  Values  Counted  Using  Equations  (12) 


^O^cp 


I  I  I 


I  ^  I  v^: 


<n 


q  cno\^  ^ 

H  Vn  UMfN-Sl* 

d  ^  CO  00  m 

o 


.1 


of 


ai  a\0 
VO  irv  m 
CO  CO  CO 


I  I  I 


oi  ca  Cl  01  €U 


irv 

ir\ 

if\ 

o 

iO  OI 

^7  ^ 

out 


oooooooooo 

Q>Q\QsQ\Q\QNO\CJ\Cr\0\  IS-J* 

I  I  I  I  I  I  I  I  I  I  I  I 

S02P5?£r>H-4-cr>o\Om»Acn 

^  Ji-  cococa  ca  cnJnca  OI  foctt 


ITkOV 

cp  Q\b  irv^ 
On  CO'S  POO 

«  •  •  •  - 

H  H  O  O 


^[CO  ifNH  O  O  O  O 


o 


.^OOOHHHHHHOIOl 
irMrvlfMf\lfMf\tfMr\lfMfMf\tf\ 


0) 


I 

Oi 


s 

H 

I 


O 

n 


It 

ii 

H  5 


.1 

O  0 

,?  JP 

Oi  Oi 


^3 


0 


-  log  A 


2  ^TCViP  5  tBlftjg 

—  log  - - —  -  —  log  M*  +  log  - - . 

3  3«  3  6 


-B 


L  V®o  J 


(40) 


Equation  (40)  «m  t«tt«d  using  eonstaats  for  Equation  (39)  previously 
dsterained  (sse  Mble  lU)  on  sBOOthed  n'^lkane  kineaatie  viscosity 
dataT. 


Values  of  Tg  'ware  calculated  using  Equation  (40)  for  a  group 
of  n-aUcanes  ranging  froe  2  equals  7  to  64  carbon  atoas  vitb  t«100  seconds 
aikfsa*^.  The  calculated  values  are  presented  in  Table  XVI.  The  values 
uere  plotted  according  to  the  aethod  of  Fox  and  Loshaek^^o^ 


1 


■•g 


(41) 


to  yield  Tg^,  defined  as  the  glass  transition  ta^peimture  at  infinite 
■oleenlar  mi^ht.  The  plot  showed  an  even'^odd  alternation  (see  Table  XVI) 
which  aade  it  difficult  to  see  trends  in  the  data.  This  effect  was 
reaoved  by  plotting  average  values  for  each  pair  of  data  points^ 


where  n  is  an  odd  Integer.  The  result  is  shown  in  Figure  6. 
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TABIE  XVI 


Tg  Values  Ccnpated  for  n-AUcanes  Using  Eguaticm  (39) 


T«-To 

67.9 

13.4 

75.0 

16.1 

88.4 

17.1 

89.4 

20.3 

102.1 

20.9 

102.9 

23.5 

1Q9.9 

24.5 

111.4 

26.5 

115.6 

27.5 

117.2 

29.0 

121.5 

29.7 

122.6 

30.9 

127.9 

31.0 

128.5 

32.1 

133.6 

32.0 

135.8 

32.7 

175* 

41.* 

♦  Based  on  data  of  A.  K.  Doolittle,  Ref  8.  Equation  (39)  paraaeters 
shown  in  Table  IH, 
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1 

Extrapolation  of  l/Tg  to  —  *  0  cannot  be  performed  accurately 

because  of  the  curvature  of  the  function  and  the  lack  of  higher  alkane 
data.  It  does  appear  that  l/Tg®^  vill  be  near  k  x  10-3,  however^  giving 
a  value  of  Tg^  of  approxinately  250*K.  This  corresponds  closely  to  the 
n  transition  range  in  polyethylene,  vhich  lies  near  250*K, 
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6,  Traaaltlon  State  TheraodynMn^  <^«  •tn  the  Glass  Transition  Region 

In  the  study  of  the  n^-aUomes  described  in  Section  3,  it  vas 
observed  that  (Tg-To)  increases  vith  Molecular  weight.  In  an  effort 
to  explain  this,  we  examined  the  relationship  between  the  Fulcher 
equation  and  the  slji5)le  Arrhenius  equation  for  the  viscosity  of  liquids. 
According  to  the  two  expressions,  the  following  equality  exists: 


log  77  *  log  A  +  B/(T-To)  **  log  A'  +  - ^ -  (43) 

2.303  BT 

Take  the  derivative  of  log  7^  with  respect  to  l/T  and  rearrange  to  obtain 

2.303  SB  1-^  ^  -I  ,  (Ut) 


the  activation  energy  in  terms  of  the  Fulcher  parameters.  At  the  glass 
tesperature  this  becomes 


E 


g 


2.303  BB 


T  “T 


(45) 


The  use  of  kinematic  viscosity  rather  than  absolute  viscosity  leads  to 
a  complication.  If  kinematic  values  are  used.  Eg  in  Equation  (45)  gives 
the  kinematic  value  of  Eg.  Differentiation  of  Equation  (39)  yields  the 
resxilt 


Eg,  absolute  Eg,  kinematic  +  2.303B!0(Tg^,  (46) 


where  Of  is  the  expansion  coefficient.  CX  is  equal  to  1  x  10“3  deg"^  to  a 
good  approximation  for  the  n-alkanes,  leading  to  a  correction  term  of 
0.04  kcal.  for  a  glass  tenperature  of  100®K,  The  calculated  values  of 
Eg  kin.  and  Eg  abs.  are  given  in  Table  3CVII. 
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Activation  Energies  for  n^-Alkanes 


NutBiber  of 

Carbon  Atoos 

Eg  kin., 
kcal. 

5 

3.59 

6 

4.01 

7 

4.80 

8 

4.92 

9 

5.65 

10 

5.75 

11 

6.17 

12 

6.26 

13 

6;  54 

14 

6.66 

15 

6.92 

16 

7.01 

17 

7.34 

18 

7.38 

19 

7.70 

20 

7.82 

64* 

10.3* 

*  Based  on  data  of  A.  K,  Doolittle, 
shown  in  Table  HI. 


Eg  abs. 


-  X  1 

kcala 

/ 

3.61 

5.32 

4.04 

5.39 

4.84 

5.48 

4.96 

5.52 

5.70 

5.58 

5.80 

5.64 

6.23 

5.67 

6.32 

5.67 

6.60 

5.-71 

6.72 

5.73 

6.99 

5.75 

7.08 

5.77 

7.41 

5.80 

7.46 

5.81 

7.78 

5.82 

7.90 

5.82 

10.4* 

5.95* 

Ref.  8.  Equation  (l4)  paraaaeters 


Bxawl nation  of  the  data  revealed  that  the  ratio  Eg  ahs./^g  Is 
nearly  constant  over  the  range  of  tdkanes  studied^  and  ajipears  to  he 
approaching  a  trtOy  constant  value  at  higher  molecular  velghts. 

This  treatment  vas  not  based  on  experimental  values  of  Tg. 

Therefore,  a  study  of  literature  data  on  the  tenperatTore  dependence  of 
viscosity  and  dielectric  relaxation  time  for  a  variety  of  polymers  vas 
Initiated,  The  bulk  of  the  viscosity  data  vas  obtained  by  Fox  and  Flory50^51. 
The  dielectric  dispersion  data  vas  published  by  Salto^3. 

It  has  been  aoply  demonstrated  In  the  literature  that  the 
temperature  dependence  of  the  shear  viscosity  and  dielectric  relaxation 
time  of  amorphous  polymers  can  be  expressed  over  broad  temperature  ranges 
above  the  transition  region  by  the  Fulcher  equation57#59>ot)^  Miller 57 
has  shown  that  the  veil-known  WIF  equation  is  mathematically  equivalent  to 
the  Fulcher  equation.  Cohen  and  Tumbull33  have  derived  a  slightly 
different  relationship  for  self-diffusion  using  a  statistical  mechanical 
approach: 


log  D  «  log  A  +  log  T^/2  + -  .  (47) 

T-To 

The  utility  of  this  equation  has  been  discussed  by  MUler^*^.  Angell^ 
has  bad  considerable  success  with  It  in  describing  the  temperature 
behavior  of  self-diffusion  and  electrolytic  conductance  in  glass-forming 
molten  salt  mixtures.  It  vas  shown  In  Section  1  that  another  relationship. 


B 

log  71  =  log  A  -  log  T  +  -  , 

T-To 


m 


fits  n-alkane  data  better  than  the  Fulcher  equation,  as  Indicated  by 
much  lower  variance  estimates  In  lecust-squasres  calculations.  Equation  (48) 
is  a  very  convenient  one  for  the  study  of  the  transition  state.  This 
will  now  be  demonstrated. 

The  Eyring^  theory  of  the  transition  state  defines  the  free  energy 
of  activation,  i.e.,  the  free  energy  difference  between  the  initial  state 
and  the  transition  state  in  the  expression 


50 


2.303  RT 


(49) 


log  ks  =  log 


/Ck 

h 


+  log  T 


> 


lAiere  ka  is  the  reaction  rate,  /C  is  a  transMssion  factor  usually 
taken  as  unity,  vhich  accounts  for  the  fact  that  not  «QJL  species  reaching 
the  transition  state  actually  decompose  to  products,  and  the  other  symbols 
have  their  usual  meanings.  This  theory,  althou^  originally  conceived 
to  describe  chemical  reactions,  has  frequently  been  applied  to  other 
irreversible  phenomena  including  viscous  flow9  and  dielectric  relaxationPl. 

The  description  of  viscosity  and  dielectric  relaxation  time 
requires  the  inversion  of  Equation  (49),  because  these  properties  are 
inversely  proportional  to  late  or  frequency  of  motion: 

log  7)  »  log  A  -  log  T  +  - - —  ,  (50) 

2.3©3  RT 


h 

Here  the  - - -  term  is  sinply  called  A,  indicating  that  it  is  a  coxmtant. 

<k 

It  is  felt  that  A  should  be  a  constant  independent  of  teapexature,  since  the 
only  variable  is  K  <U3d  this  should  depend  only  on  the  nature  of  the 
polymer. 

Coaparison  of  Equation  (48)  and  (50)  indicate  that  as  T  approaches 
infinity,  the  two  equations  becoam  identical,  with  A®^ /2»303  R  «  B.  This 
means  that  the  values  of  A  in  both  equations  are  the  same,  leading  to  the 
following  result  for  all  teaperatures: 


Ag+ 


2.303  RBT 

T-To 


(51) 


The  activation  energy  (or  E/^)  is  obtained  by  differentiation  of  log  7^ 

with  respect  to  — ,  giving 
T 


RT  +  2.303  RB 


2 


T-To  I 


(52) 


51 


enthalpy 


.Finally^  liy  aalrlTiB  the  xuraal  assvm^rtion  that  the  activation 


is  related  to  Ey|  hy 


-  RT, 


(53) 


ve  can  obtain  the  entropy  of  activation: 


T 


m 


Equations  (51)  and  (52)  shew  that  the  activation  energy  and 
free  energy  increase  with  decreasing  tenperature,  approaching  infinity 
as  T  approaches  To.  It  can  he  shewn  that  the  same  is  true  of  the 
entropy  of  activation.  Ccoibination  of  Equations  (5l)^  (52)^  and  (54) 
yields  the  expression^ 


^  2.303  KB 

- 

T-To 


(55) 


idiich  clearly  goes  to  infinity  at  Tq. 

The  viscosity  and  relaxation  time  data*^3>50-52  used  to 
evaluate  the  parameters  A,  B,  and  Tq  in  Equation  (48)  enploying  a 
non-linear  least-squares  method.  The  transition  state  functions  E^(and 
E/^),  and  were  then  conputed  at  the  glass  temperature.  ‘ 

Tables  XVm  and  XIX  give  the  results  of  these  calculations. 

Eyring9  derived  an  alternate  expression  for  the  free  energy  of 
activation  in  viscous  flow  by  assuming  that  the  molecules  of  a  liquid 
move  from  one  equilibrium  position  to  another  under  the  influence  of  an 
applied  force: 


AGv  ■  2.303  RT  [log  7)  +  log  V  -  log  Kh].  (56) 
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TABIE  mil 

Transition  State  Parameters  at  Tg  for  Dielectric  Relaxation,  Based  on  Saito*s^3  Data 
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*  Tg  calculated  using  Equation  and  Tg  values  of  Saito, 
**  Variance  estimate,  least  squares  fit. 


The  subscript,  v,  indicates  that  the  fJree  energy  defined  by  Equation  (56) 
is  being  used.  By  con5)aring  activation  enthalpies  of  linear  polyesters 
irith  heats  of  vaporization  of  analogous  sis^le  esters,  Syriug  concluded 
that  the  molar  volume  V  of  a  moving  polymer  segment  is  about.  500  cbuS.  In 
order  to  follow  his  theory  throu^,  we  shall  assume  the  same  value. 

Equations  (56)  and  (^)  may  be  conibined  to  yield 


Ag^  *  2.303  RT  [log  A  -  log  T  + - -  +  log  V  >  log  Mh].  (57) 

T-Tq 


Now  AGt> 

of  the  brackets  to  give 


A®v  "  AOt 


Equation  (58) 
of  log  V  and  log  Nh, 


Ag^  -  A(4 


free  energy  defined  in  Equation  (5I)#  is  talcen  out 

2.303  RT  [log  A  -  log  T  +  log  V  -  log  Nh].  (58) 

is  simplified  numerically  by  substituting  the  value 

2.303  RT  [log  A  -  log  T  +  5.098].  (59) 


Values  of  A^^  a°d  As^  computed  using  Equation  (59)  are  shown 
in  Table  XX.  It  me^  be  seen  that  the  free  energy  values  are  somewhat 
higher,  and  the  entropy  values  slightly  lover  than  those  in  Table  XIX. 
Thus,  the  results  are  still  basically  the  same  as  before.  Therefore,  the 
discussion  which  follows  applies  to  both  sets  of  results. 


It  is  quite  obvious  that  in  contrast  with  the  n-alkane  results 
in  Table  XVII  the  ratio  of  activation  energy  to  tei^erature  is  not  a 
constant  at  the  glass  temperature,  except  possibly  within  the  same  polymer 
series. 


Several 
glass  temperature 


comments  on 
can  be  made. 


the  nature  of  the  transition  state  near  the 
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lABIB  XX 

Transition  State  Paraoeters  at  Tg  Based  on  Eyrlng  Viscosity  Theory 


Polyaer 

Mria  X  l<r^ 

Icc&X* 

Sy,  «a 

PS 

1.930 

31.5 

564 

1.330 

30.5 

533 

0.665 

32.5 

74i 

0.359 

26.2 

451 

PIB 

0.54 

17.2 

143 

Becaaethylene  Adipate 

0.451 

16.2 

129 

Diethylene  Adipate 

0.P90 

17.6 

481 
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First,  the  activation  energy  is  quite  large,  around  200  kcal., 
at  the  glass  ten5>erat\jre,  except  for  those  polymers  with  very  low  glass 
temperatures.  It  is  obvious  from  this  that  a  rather  large  region  in 
the  polymer  must  be  involved  in  motions  leading  to  the  glass -transition. 

The  activation  entropy  is  large  and  positive  in  every  case. 

This  indicates  that  the  formation  of  the  transition  state  requires 
considerable  disordering  of  the  polymer  in  the  transition  region.  At 
least  two  explanations  for  this  are  possible.  First,  the  size  of  the 
volume  element  affected  by  the  motion  of  a  particular  group  may  increase, 
and  second,  there  may  be  an  increase  in  short-range  order  in  the  neighbor¬ 
hood  of  the  transition  temperature  necessitating  the  creation  of  relatively 
more  disorder  to  reach  the  transition  state.  It  is  likely  that  both  effects 
contribute  to  the  large  positive  entropy  of  activation. 

The  free  energy  of  activation  at  the  glass  temperature  was  found 
to  be  positive  in  every  case.  The  variation  in  this  property  was  much  less 
than  that  of  the  entropy  or  activation  energy,  ranging  from  a  value  of  12 
for  poly(decamethylene  adipate)  to  30  for  PMMA, 

A  rather  obvious  result  of  this  work  is  the  close  cozrespondence 
between  the  (X  dielectric  dispersion  process  and  the  viscous  flow  process 
above  Tg.  Althou^  both  processes  were  not  studied  for  any  given  polymer, 
the  activation  energies  seem  to  be  approximately  equal  for  polymers  with 
the  same  Tg  in  either  process. 

This  correspondence  between  dielectric  and  mechanical  properties 
has  been  observed  before^S,  TheSe  results  supply  additional  evidence  for 
tbe  equivalence  of  dielectric  and  mechanical  relaxation  phenomena.  Althou^^ 
the  similarity  between  the  response  of  a  dipole  to  an  alternating  electric 
field  and  the  xesponse  of  a  polymer  to  a  mechanical  stress  may  not  be 
appeu^nt  at  first,  it  is  easily  understood  when  we  consider  the  size  of 
the  region  involved  in  the  response  of  a  single  dipole.  The  magnitude  of 
the  activation  energy  at  Tg  is  a  clear  indication  that  this  region  is 
quite  large,  involving  coordination  with  seemingly  unrelated  portions  of 
the  polymer.  This  is  qualitatively  similar  to  the  coordination  required 
for  Hiechanlcal  relaxation. 

It  has  been  observed  by  several  workers^3^64  ^that  the  activation 
energy  for  the  (X  relaxation  process  has  a  maximum  value  near  the 
dllatometric  glass  transition  temperature.  Figure  7  shows  the  results 
obtained  by  Sommer®3  from  dynamic  mechanical  measurements  on  poly( vinyl 
chloride)  with  an  average  molecular  weight  of  35,000.  Superimposed  is  the 
activation  energy  curve  predicted  by  Equation  (48)  and  (52),  based  on 
Saito*s^3  dielectric  data  on  poly(vii^l  chloride)  with  equal  to  36,000. 
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It  can  be  seen  that  the  yalne  of  the  activation  energy  at  Tg  is  someirhat 
higher  for  the  mechanical  process  than  the  dielectric  process.  This  is 
not  \mreasonable,  especially  in  the  case  of  polymers,  idiere  the  coordination 
required  for  motion  of  a  single  dipole  should  be  more  localized  than  that 
required  for  over-all  mechanical  relaxation.  If  this  rationale  is  correct, 
the  dielectric  and  mechanical  activation  energies  should  be  nearly 
identical  for  polar  molecxiles.  Davidson  and  Cole°®  found  this  to  be 

so  in  the  case  of  n-propanol  and  glycerine.  Actually  these  authors 
evaluated  the  Fulcher  equation  constants  for  dielectric  relaxation  time 
and  viscosity  and  found  B  and  Tq  to  be  virtually  identical  for  the  two 
processes,  which  of  cotirse  indicates  equal  activation  energies. 

Below  the  glass  tai^erature  the  experimental  (mechanical)  value 
of  the  activation  energy  declines,  but  the  calculated  value  continues  to 
increase  and  will  in  fact  go  to  infinity  at  Tq. 

How  can  the  discrepancy  between  experimental  observation  and 
calc\ilation  be  explained?  Fox  and  Floryl5  suggested  that  the  decirease  in 
activation  energy  was  due  to  a  condition  of  constant  free  volume  below  Tg. 

A  more  satisfactory  eaplanation  has  been  proposed  by  Saito^3.  Hie  found 
that  the  activation  energy  for  dielectric  relaxation  can  be  increased  in 
the  transition  region  by  cooling  the  sample  to  the  measurement  temperature 
more  slowly.  He  proposed  that  the  decrease  in  activation  energy  is  caused 
by  the  non-equilibrium  nattire  of  the  sample  below  Tg.  In  other  words,  the 
time  required  for  equilibration  below  Tg  is  so  long  that  the  polymer 
behaves  eui  though  it  were  still  at  some  higher  temperature. 

The  activation  energy  peak  in  Sommer's  plot  occurs  at  3^1 ’K, 
which  is  within  the  range  reported  in  the  literature  for  the  Tg  of 
poIy( vinyl  chloride).  It  would  appear  on  the  basis  of  the  discussion 
above  that  the  position  of  the  activation  energy  maxiimim  corresponds 
closely  to  the  glass  tesperature  determined  dilatcaaetrically.  The  glass 
temperature  xoeasured  in  this  way  depends  only  on  the  sasple  and  not  on  an 
arbitrary  frequency  of  measurement. 

An  illustration  of  the  proposed  method  is  shown  in  Figure  8.  The 
glass  temperature  is  shown  in  the  figure  as  the  first  tempereture  at  which 
the  e3q>erimental  data  deviates  from  Equation  (^),  since  this  is  the  first 
indication  of  a  non-equllibritim  condition  in  the  polymer.  A  hypothetical 
maximum-loss  frequency  curve  is  also  shown  in  Figure  8,  to  show  that  it  is 
not  absolutely  necessary  to  calculate  and  plot  the  activation  energy  in 
order  to  determine  Tg. 

It  would  be  of  interest  to  study  the  y5  transitions  of  several 
polymers  in  order  to  see  if  there  is  a  corresponding  peak  in  the  activation 
energy.  If  such  a  peak  exists,  it  should  lead  to  some  new  thoughts  con¬ 
cerning  the  nature  of  the  tiransition,  particularly  regarding  free  volume. 
Similarly  it  should  be  detonnined  if  there  is  a  secoM  in  the 

activation  energy  of  the  fj  process  in  the  neighborhood  of  Tg. 
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VII.  IntermolecTilar  Interactions 

The  forces  of  interaction  between  molecules  play  an  important 
role  in  the  equilibrium  and  transport  properties  of  liquids,  and  hence 
also  in  the  transition  to  other  states  such  as  the  glass.  It  is  not 
STirprising,  therefore,  to  find  that  many  theories  and  enqpirical  cor¬ 
relations  related  to  the  glass  transition  temperature  are  expressed  in 
terms  of  some  measure  of  intermolecular  attraction,  such  as  the  cohe¬ 
sive  energy  density,  the  hole  energy  or  the  thermal  coefficient  of 
expansion  of  free  volume .  The  press^lre  coefficients  in  Table  VI  are 
also  a  measure  of  intermolecular  forces .  It  would  be  desirable  to  be 
able  to  predict  the  effects  of  molecular  stmicture  on  intermolecular 
interactions . 

A  correlation  proposed  by  Smll^^,  derived  from  the 
van  der  Waals  equation,  treats  (EV)^/^  as  a  constitutive  property, 

E  being  the  cohesion  energy  per  mole  and  V  the  molar  volume.  Since 
E  is  a  function  of  temperature,  an  attempt  was  made  to  adapt  this 
correlation  to  the  prediction  of  the  Antoine  parameter  B  ,  which  is 


temperature independent,  by  means  of  the  relation 


E  =  2.303R  BpT^/(T-Tq)^-RT. 


(60) 


The  values  of  the  Antoine  constants  Bp  and  T^j  selected  by  API  Project 
for  the  n-alkanes'  were  submitted  to  the  Mathematical  Analysis 
Group  to  test  the  function 


”1 


2.303  EBpT^/(T-To)  -  m  V^/^  = 


ai  n  +  bi 


.1/2 


which  follows  frcm  Equation  (60)  and  the  additivity  of  (EV)  '  .  The 

values  of  the  adjustable  parameters  a^^  and  b^  were  chosen  so  as  to 
minimize  the  sum  of  squares  of  the  deviations  in  B  .  As  shown  in 
Table  XXII,  the  fit  to  the  data  is  quite  poor,  which  may  be  due  to  the 
approximations  inherent  in  Equation  (60),  rather  than  shortcomings 
in  Small’s  correlation. 

1/2 

It  was  thoxight  that  (B  V)  might  be  an  additive  property, 
by  analogy,  but  this  is  not  the  Case.  However,  it  has  been  found  that 


(BpV) 


0.7240 


=  835-59  n  +  732.14 


(62) 


gives  a  very  good  fit,  as  seen  in  Table  XXIII .  A  number  of  empirical 
correlations  for  the  Antoine  constants  were  stxidied,  and  the  results 
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TableXXI 


Constants  for  the  n-Para££in  Hydrocarbons* 
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!  4 
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945.90 
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0.62143 
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1064.63 
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0.65486 

48.794 

1171.530 

1 

147 .456 

0.67957 

56.260 

1268.115 

8 

163.530 

0.69854 

63.643 

1355.126 

*  0 

179.670 

0.71386 

71.541 

1428.811 

195.905 

0.72631 

78.680 

1501.268 

} 

ili 

212.217 

0.73658 

85.138 

1572.477 

|i2 

228.579 

0.74522 

92.849 

1625.928 

p 

244.924 

0.75276 

98.877 

1689.093 

\.lr 

261.312 

0.75923 

105.626 

1739.623 

s 

277.698 

0.76494 

111.869 

1789.658 

1 

;io 

294.083 

0.77002 

118.632 

1831.317 

L 

/ 

310.51 

0.77445 

123.140 

1880.61 

L8 

326.93 

0.77847 

128.63 

1920.60 

19 

343.36 

0.78207 

133.56 

1961.6 

20 

359.83 

0.78525 

139.96 

1994.0 

Values  taken  from  Reference  T 


Table  XXII 


\  -  [s  +  bi)"  +  592.48]j<2-2|j|^)'] 

=  1.5845  X  10^ 

=  1.7165  X  10 


n 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


Calculated  B 

- p 

605.24 

823.53 

1012.43 

1178.736 

1322.793 

1444.803 

1538.378 

1620.845 

1694.582 

1730.077 

1777.412 

1796.684 

1811.130 

1802.062 

1826.68 

1820.97 

1817.6 

L772.1 


(Observed  B^-  Calculated  B^) 
207.96 
122.37 
52,20 
-7.206 
-54.678 
-89.677 
-109.567 
-119.577 
-122.105 
-104.149 
-88.319 
-57.061 
-21.472 
29.255 
53.93 
99.63 
144.0 
221.9 
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Table  XXIII 


_  dr  u 

“  m  ^^2"  +  ^2^ 

Eg  -  8.3559  X  10^ 
bg  »  7.3214  X  10^ 
e  g  •=  1.3813 

n  Calculated  B  (Observed  B  -  Calculated  B  ) 

—  - p  .i - =c.ir - p - p^ 


3 

789.03 

24.17 

4 

955.72 

-9.82 

5 

1080.32 

rl5.69 

6 

1184.205 

-12.675 

7 

1274.868 

-6.753 

8 

1355.915 

-0.789 

9 

1430.191 

-1.380 

10 

1498.587 

2.681 

11 

1562.112 

10.365 

12 

1621.634 

4.294 

13 

1678.130 

10.963 

14 

1731.577 

8.046 

15 

1782.607 

7.051 

16 

1831.474 

-0.157 

17 

1878.14 

2.47 

18 

1923.09 

-2.49 

19 

1966.4 

-4.8 

20 

2008.0 

-14.0 

6k 


are  shown  in  Tables  XXIV  to  XXVII.  The  three  predictors  for  Tq  are  roughly 
equivalent,  and  that  shown  in  Table  XXVII  is  probably  the  best  choice 

(To/d)  =  146.879  n  +  27.51I+  (63) 

jEquatlon  (62)  shows  that  a  constant  increment  of  835*59  in 
(B-nV)  *  ’  results  from  an  increase  in  chain  length  of  one  atom.  Re- 
arranging,  and  defining  a  new  parameter  Q  =  (BpV)  ‘  /835.59#  have 

Q  =  n  +  0.876  (61^) 


for  n-alkanes.  Q  was  calcxjuLated  for  some  branched  alkanes,  using 
values  of  Bp  and  V  tabiilated  by  API  Project  It  was  found  that 

the  presence  of  a  branch  reduced  the  valtie  of  Q  by  about  0.15  for 
each  side  group  s,  compared  with  the  e35)ected  value  for  the  normal 
isomer  of  the  same  n.  To  illustrate  the  magnitude  of  the  deviations 
from  this  rule,  the  values  of  Q  -  n  +  0.155  "  O.876  are  listed  in 

Table  XXVUI. 

The  valxjes  for  2,3-dimethyl  isomers  are  omitted  from 
Table  XXVIII  because  they  were  higher  by  an  average  of  0.037  than  those 
calculated  from  Q  -  n  +  0.15s  -  O.876.  While  the  deviation  is  very 
small,  it  seemed  to  be  consistent,  so  a  number  of  other  compounds 
with  vicinal  side  chains  were  studied.  The  average  deviation  for 
12  compounds  with  vie  side  chains,  and  for  which  Bp  was  reported  to 
at  least  0.1*,  was  0.035  for  each  pair  of  vie  side  groups  v.  To 
illustrate  the  adherence  to  this  rule,  values  of  Q  -  n  +  0.15j^  -  0.035  v 
-  0.876  are  listed  in  Table  XXIX^  Data  for  the  last  eight  compounds 
in  the  table  are  from  Driesbach.^ 

The  largest  errors  in  the  use  of  Q  =  n  -  0.15s  +  0.035  ▼ 

+  0.876  occur  in  those  compounds  where  Bp  is  reported  only  to  the 
nearest  degree,  and  are  therefore  the  least  reliable.  The  worst  case 
is  3 -isopropyl-2 -methylhexane.  Table  XXIX,  for  which  the  calcuJ^ted 
Bp  =  IU7U®,  an  error  of  45*  or  3^^-  Of  those  compounds  for  which 
Bp  is  reported  to  better  than  0.1*,  the  largest  error  occurs  with 
2-methylheptane,  for  which  the  calculated  Bp  »  1315*96*,  an  error  of 
21.51*  or  1.65^.  The  most  convincing  evidence,  however,  of  the  pre¬ 
cision  of  this  correlation,  is  the  fact  that  it  is  able  to  distin- 
quish  a  third-order  effect,  that  due  to  the  placement  of  side  chains 
on  adjacent  carbon  atoms. 

It  is  reemnmended  that  these  correlations  for  the  Antoine 
constants  of  alkanes  be  extended  to  include  other  classes  of  coin5>ounds. 
Furthermore,  the  structural  effects  which  have  been  observed  must  have 
an  Influence  on  glass  formation,  although  this  aspect  of  the  problem 
has  not  yet  been  adequately  studied. 
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Table  XXIV 


n 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


Bp  =  d[»^n  +  b3] 

a^  ■=  8.1227  X  10^ 
b3  =  3.7871  X  10^ 

-  4.1048  X  lO”^ 


Calculated  B 

- p 

779.74 
953.52 
1081.44 
1186.632 
1277.565 
1358.330 
1432.113 
1499.958 
1562.925 
1621.945 
1678.023 
1731.141 
1781.934 
1830.668 
1877 .27 
1922.30 
1965.7 
2007.5 


(Observed  B^-  Calculated  B^) 
33.46 
-7.62 
-16.81 
-15.102 
-9.450 
-3.204 
-3.302 
1.310 
9.552 
3.983 
11.070 
8.482 
7.724 
0.649 
3.34 
-1.70 
-4.1 
-13.5 
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Table  XXV 

To  =B(a’i+b)''4 
P  H 

a,  =  7.6210  X  lO”^ 

4 

b,  =  3.7362  X  lO"^ 

4 

c,  =  4.07977  X  lo"^ 

4 

Calculated  (Observed  T/>  -  Calculated  \ 


3 

24.73 

0.43 

4 

32.98 

0.18 

5 

41.12 

0.04 

6 

49.099 

-0.305 

7 

56.889 

-0.629 

8 

64.442 

-0.799 

9 

71.501 

0.040 

10 

78.611 

0.069 

11 

85.768 

-0.630 

12 

92.033 

0.816 

13 

98.914 

-0.037 

14 

105.120 

0.506 

15 

111.341 

0.528 

16 

117.073 

1.559 

17 

123.328 

-0.188 

18 

129.01 

-0.38 

19 

134.78 

-1.22 

20 

139.98 

-0.02 

6? 


Table  XXVI 


T 

o 


=  tv 


5 


a^  -  2.3518  X  10 
=  7.8325 
-  1.7555 


n  Calculated  T«  (Observed  T^  -  Calculated 


3 

23.64 

1.52 

4 

33.05 

0.11 

5 

41.58 

-0.42 

6 

49.598 

-0.804 

7 

57.257 

-0.997 

8 

64.619 

-0.976 

9 

71.760 

-0.219 

10 

78.689 

-0.009 

11 

85.430 

-0.292 

12 

92.007 

0.842 

13 

98.461 

0.416 

14 

104.779 

0.847 

15 

110.989 

0.880 

16 

117.099 

1.533 

17 

123.098 

0.042 

18 

129.014 

-0.384 

19 

134.844 

-1.284 

20 

140.583 

-0.623 
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Table  XXVII 


To  “  +  b^J  ^ 

a,  =  4,6879  x  10 
6 

=  2.7514  X  10 

c.  =  7,5484  X  lO"^ 
o 


n  Calculated  Tq  (Observed  Tq  -  Calculated  T/,  ) 


3 

23.59 

1.57 

4 

33.03 

0.13 

5 

41.56 

-0.40 

6 

49.596 

-0.802 

7 

57.262 

-1.002 

8 

64.627 

-0.984 

9 

71.771 

-0.230 

10 

78.702 

-0.022 

11 

85.442 

-0.304 

12 

92.018 

0.831 

13 

98.471 

0.406 

14 

104.787 

0.839 

15 

110.994 

0.875 

16 

117 . 100 

1.532 

17 

123.095 

0.045 

18 

129.01 

-0.38 

19 

134.83 

-1.27 

20 

140.57 

-0.61 
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Correlation  for  Intennolectilar  Interactions 


Effect  of  Vlcliaal  Side  Groups  on 
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